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We determine the quark masses and the chiral condensate in the MS scheme at NNLO from Lattice QCD in 
the quenched approximation at f3 = 6.0, f3 = 6.2 and /3 = 6.4 using both the Wilson and the tree-level improved 
SW-Clover fermion action. We extract these quantities using the Vector and the Axial Ward Identities and non- 
perturbative values of the renormalization constants. We compare the results obtained with the two methods and 
we study the 0(a) dependence of the quark masses for both actions. 



1. INTRODUCTION 

To describe the light hadron spectrum one con- 
siders massless QCD a good approximation to 
this world and the vacuum not symmetric un- 
der chiral transformations. The chiral condensate 
is the order parameter which governs the spon- 
taneous chiral symmetry breaking. The pseu- 
doscalar mesons are identified with the goldstone 
bosons and their physical masses are attributed 
to the mass term in the QCD Hamiltonian. Since 
free quarks are not physical states, quark masses 
cannot be measured directly in the experiments 
and can be determined from the meson spectrum 
using non-perturbative techniques. On the lattice 
one can compute the quark masses and the chiral 
condensate from first principles and it is the only 
procedure that can be systematically improved. 
The methods and the symbols we have used and 
all the results we have obtained are fully described 
in II- 

2. QUARK MASSES 

The usual on-shell mass definition cannot be 
used for quark masses and their values depend on 
the theoretical definition adopted. In the follow- 
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ing we will give our final results for the running 
quark masses defined in the MS scheme. Quark 
masses can be defined from the Vector Ward Iden- 
tity (VWI). Neglecting terms of 0(a), the VWI 
between on-shell hadronic states can be written 
as [| 

(*\r%W = l(±-±)(a\sw>. (i) 

Eq. (Q) fixes the relation between the lattice bare 
quark mass in lattice units and the hopping pa- 
rameter, i.e. for the Wilson action m = l/2(l/fc— 
l/k c ). Quark masses can also be extracted from 
the Axial Ward Identity (AWI). Neglecting terms 
of 0(a), the AWI can be written as ||] 

Z A (a\d^A a ^) = 2(m ~m){a\P a \P) , (2) 

where m is defined in || . The light and strange 
quark masses are determined by fixing to their 
experimental values the masses of the tt and K 
mesons. The standard perturbative approach 
|l],^| uses lattice and the continuum perturba- 
tion theory to connect the "bare" lattice quark 
mass to the renormalized MS one. The scale 
1/a, where a is the lattice spacing, of our simula- 
tions is a -1 ~ 2 — 4 GeV. At these scales we ex- 
pect small non-perturbative effects. However the 
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Table 1 

Summary of the parameters of the runs analyzed in this work. 

Matrix Elements 

C60a C60b W60 C62 W62 C64 W64 

P 6.0 6.0 6.0 6.2 6.2 6.4 6.4 

Action SW SW Wil SW Wil SW Wil 

# Confs 490 600 320 250 250 400 400 
Volume 18 3 x 64 24 3 x 40 18 3 x 64 24 3 x 64 24 3 x 64 24 3 x 64 24 3 x 64 

Renormalization Constants 

C60Z W60Z C62Z W62Z C64Z W64Z 

P K0 K0 6^2 6^2 6A 6A 

Action SW Wil SW Wil SW Wil 

# Confs 100 100 180 100 60 60 
Volume 16 3 x 32 16 3 x 32 16 3 x 32 16 3 x 32 24 3 x 32 24 3 x 32 



Table 2 

Quark Masses an the Chiral Condensate from the 
VWI in MeV. WS values are at p = 2 GeV. 



Run 








W60 
W62 
W64 


5.8(2)(1) 
5.4(2)(1) 
4.9(2)(1) 


130(3)(2) 
124(4) (2) 
112(5)(2) 


(247±2± l) 3 
(250±3±1) 3 
(258±4± l) 3 



"tadpole" diagrams, which are present in lattice 
perturbation theory, can give rise to large pertur- 
bative corrections and hence to large uncertain- 
ties in the 1-loop matching procedure at values 
of P — &/g\ — 6.0 — 6.4. A non-perturbative 
renormalization (NP) technique eliminate these 
uncertainties |5|||. Za, Zp and Z$ can be cal- 
culated by imposing the renormalization condi- 
tions, proposed in on the quark states of mo- 
mentum p 2 = p 2 and in the Landau gauge gj. 
This procedure works if p satisfies the condition 
Aqcd •C p <C 1/a to avoid chiral symmetry 
breaking effects, large higher-order perturbative 
corrections and discretization errors. In figure |l| 
we show the scalar renormalization constant di- 
vided by its Renormalization Group (RG) evolu- 
tion (Zg GI ). The data show that discretization 
errors are within statistical errors in the range 
0.5 < pa < 2 where the NP renormalization is 
applied to compute the quark masses. The NP 
method allows a fully non perturbative defini- 
tion of the renormalized quark masses in the RI 



scheme 



i ri {p) = 



ma 



m H1 (//) 



yRI 

^ A — 1 

A -pa 1 



where for large time separations 

, s 1 • , n(^o(t)P(0)) 
p(a) = - sinh(Mpg) V ; / ; // 



(3) 



(4) 



The MS definition of the renormalized quark 
masses is intrinsically perturbative and can be 
related to the RI one through continuum per- 
turbation theory only: 



in 



MS M = U™ s (p,p>) 



where U^ s (p, p') is the RG evolution of the 
quark mass. Z^ s /Z^ is the matching factor 
computed in perturbation theory at scales p ~ 
2 — 4 GeV large enough to avoid non-perturbative 
effects and/or higher order corrections. 

3. THE CHIRAL CONDENSATE 

The chiral condensate can be defined using the 
AWI arising from the variation of the non-singlet 
pseudoscalar density. The integrated AWI be- 
comes m 
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lim 2(r7i -m) / d 4 x(P(x)P(0)) . 
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Table 3 

Quark Masses and the Chiral Condensate from 
the AWI in MeV. WS values are at \i = 2 GeV. 



Run 



C60a 

C60b 

W60 

C62 

W62 

C64 

W64 



N 



6-0(2) (9) 
5.7(2)(8) 
5.7(2)(8) 
5.8(5)(6) 
5.4(2)(5) 
4.4(3)(3) 
4.7(2)(3) 



136(4) (20) 
132(4)(19) 
127(4) (17) 
131(7) (14) 
122(4)(12) 
104(5) (6) 
108(5)(8) 



(242±3± 12) 3 
(244±2±12) 3 
(248±2±11) 3 
(245 ±4 ±9) 3 
(251 ±3 ±8) 3 
(265 ±4 ±5) 3 
(262 ± 4 ± 6) 3 



Figure 1. Z§ GI vs fi for the run W64. 
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In the chiral limit one can show that the above 
expression is equivalent to 



PpMI 



-^-(ipip) = - lim 

JS f m -*m c 4(r7lQ — m) 



(6) 



which is the familiar Gell-Mann-Oakes-Renner 
(GMOR) relation. We write the relation (g) for 
the renormalized condensate as 



HS 



1 



1 



1 t 2 qAWI 



-a~ L f 

Nf^ 1 -" 2 J *Z A 



(7) 



where f x = 0.1282 GeV is the "experimental" 
value in physical units H and 



Mp = 



2ap 



C HS 
1 

QAWI 



1 1 



(8) 



The main advantages of the above formulas are 
to avoid the error amplification of the standard 
method because we are left with only one power 
of the UV cutoff aT 1 and to determine the slope 
(jHS piffjj without extrapolation to the chiral 
limit. Note that since we work at k values typ- 
ical of the strange quark mass, we are implicitly 
assuming that the slope will not change in the 
chiral region. 

4. RESULTS 

We have computed the renormalization con- 
stants of bilinear quark operators with the NP 



method proposed in ||. The /i-dependence of 
Za and Z$ are in excellent agreement with the 
RG predictions. Zp is in good agreement but the 
chiral symmetry breaking effects deserve a more 
accurate analysis. The parameters and the ac- 
tion used in each simulation are listed in Table 
1 and the main results we have obtained are re- 
ported in Tables | and § The data at p = 6.4 
have to be consider for an exploratory study only, 
since the physical volume and the time extension 
of the lattice are too small to be considered re- 
liable. The MS values of the quark masses and 
the chiral condensate reported in Tables § and |^ 
show a very good agreement between the values 
extracted from the VWI and from the AWI using 
the non perturbative determinations of the renor- 
malization constants, while the same comparison 
using the perturbative values of Za-, Zp and Z$ 
fails giving inconsistent results for the two meth- 
ods [Q. This pattern is also found by for 
different fcrmion actions. In the /3 range studied, 
there is no statistical evidence for an "a" depen- 
dence of the quark masses and the chiral conden- 
sate. We believe that the best estimates for the 
light and strange quark masses and the chiral con- 
densate are mf 75 ^ GeV) = (5.7±5±8±8) MeV, 
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(2 GeV) = (130 ± 8 ± 15 ± 15) MeV and 
(#) M (2 GeV) = (245±4±9±7iWW) 3 , 



where the first error is statistical, the second is 
due to the non-perturbative renormalization and 
the third is an estimate of the overall systematic 
errors on the bare quantities Q . 
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